Efficiency at maximum power of thermally coupled heat engines 
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We study the efficiency at maximum power of two coupled heat engines, using thermoelectric 
generators (TEGs) as engines. Assuming that the heat and electric charge fluxes in the TEGs 
are strongly coupled, we simulate numerically the dependence of the behavior of the global system 
on the electrical load resistance of each generator in order to obtain the working condition that 
permits maximization of the output power. It turns out that this condition is not unique. We 
derive a simple analytic expression giving the relation between the electrical load resistance of each 
generator permitting output power maximization. We then focuse on the efficiency at maximum 
power (EMP) of the whole system to demonstrate that the Curzon-Ahlborn efficiency may not always 
be recovered: the EMP varies with the specific working conditions of each generator but remains in 
the range predicted by irreversible thermodynamics theory. We finally discuss our results in light 
of non-ideal Carnot engine behavior. 

PACS numbers: 05.70.Ln, 84.60.Rb 



I. INTRODUCTION 

One of the pillars of thermodynamics is the existence 
of a limitation in the conversion of heat into work: the 
second law, as expressed by Carnot, gives the maximum 
efficiency that an ideal heat engine working between two 
temperature reservoirs Thot and Tcoid (with Thot > ^coid) 
may boast. This maximum, known as the Carnot effi- 
ciency, here denoted rjc^ is only function of these tem- 
peratures: r]c = {Thot - Tcoid)/Thot H- This limit may 
be reached only if the process of conversion is fully re- 
versible; but the ideal Carnot engine is a zero-power 
engine since a reversible transformation in a thermody- 
namic system is quasi- static: it requires an infinite time 
to complete. We can even go further by asserting that 
the Carnot engine is, as a matter of fact, unphysical: no 
dissipative element ensures its causality, so no arrow of 
time can be defined for such system. 

^From a theoretical as well as a practical viewpoint, 
only models of irreversible heat engines must be consid- 
ered to address the central points of causality and system 
evolution. The main consequence of the introduction of 
irreversibilities in models is the dependence of the en- 
gine's efficiency, r^, and produced power, P, on the work- 
ing conditions. This dependence ties 77 and P together: 
77 = r]{P). The Carnot efficiency can then be reached 
only when the engine works infinitely slowly; conversely, 
as the speed of the thermodynamic cycle of the engine 
increases, its conversion efficiency decreases. For very 
fast cycles the efficiency even reduces to zero. As regards 
power, if the cycle time r is very long, it vanishes as it is 
given by the ratio of work, finite in that case, to the time 



r. On the contrary, if the speed increases too much, ir- 
reversibilities become preponderant and power vanishes 
again. Hence, between these two extremal cases, there 
are a variety of working conditions also leading to power 
maximization. For practical purposes, real thermody- 
namic engines must produce power, not just work: this 
observation shifted interest from efficiency, at the center 
of Carnot's derivations, to output power, thus establish- 
ing the theory of finite-time thermodynamics [2|-Q ; for a 
recent review, see Ref. \§\. 

This theory, which emerged during the 1970s, may be 
viewed as a generalization of the thermodynamics of ir- 
reversible processes. The seeds were sown in the 1950s 
when physicists and engineers were working to obtain the 
best possible performances from the newly built atomic 
power plants: their design not only had to ensure power 
maximization but also an efficiency as high as possible 
for this particular working condition. To reach this pur- 
pose, a new model of heat engine was put forward: the 
irreversibilities were introduced by considering finite lin- 
ear thermal conductances between heat reservoirs and 
the ideal heat engine, i.e. with no internal dissipations. 
This model was later referred to as endoreversible system 
for irreversibilities occur outside the engine ^]. Basing 
his analysis of efficiency at maximum power [r]{Pma^)] 
on this then novel model, Novikov [7, 8] was the first to 
demonstrate the following seminal expression for EMP: 
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for endoreversible engines. Note that Novikov [7] at- 
tributes the above expression to Yvon [9] (however he 
points out that no demonstration was given). Nowadays, 
this result, very popular owing to its simplicity, is known 
as the Curzon-Ahlborn (CA) efficiency t^ca, after the 



names of the two authors who rederived Eq. ([T]) inde- 
pendently from a refined analysis of a model heat engine 
in 1975 [10|- It is interesting to note that this expression 
was derived in 1957 by Chambadal [11], who thus shares 
the paternity of Eq. ([T]); however the hypothesis made 
on the heat exchange between the heat reservoirs and 
the engine are completely different from the ones made 
by Novikov, and Curzon and Ahlborn who assumed the 
linearity of thermal contacts. Hence, even if the final re- 
sult is the same, one should keep in mind that the models 
are different. The CA efficiency turned out to be in good 
agreement with real EMP values, particularly in the case 
of power plants [10]. It has then rapidly been consid- 
ered as a universal upper bound for EMP of heat engines 

In 2005 Van den Broeck stated that the CA efficiency 
derives from the theory of linear irreversible thermody- 
namics, and that it can be generalized to any heat engine 
working in the regime of strong flux coupling [14]; one 
of the main interests of this article is the elements that 
it brings to the debates in finite-time thermodynamics. 
First, Van den Broeck used the Onsager's reciprocal rela- 
tions of the force-flux formalism p^5l-[l7| to study heat en- 
gines. Finite-time thermodynamics analyses thus moved 
from cyclic to steady-state engine operations. The cycle 
time is no longer a parameter to optimize in that case: 
power maximization is obtained by tuning the thermo- 
dynamic forces applied to the engine. One should note 
that the EMP and all features associated to heat engines 
with cyclic operations remain the same for steady-state 
operations. Second, Van den Broeck highlighted the fact 
that the CA efficiency can only be obtained for engines 
operating in the strong coupling regime. This assumption 
implies that inside the heat engine, the heat flux and the 
matter flux, which yield the output power, are tightly 
linked; more precisely they are proportional in the frame 
of linear theory. 

The use of Onsager coefficients allows a natural and 
transparent description of the internal dissipations and 
thus puts forward the question of the influence of irre- 
versibilities on EMP. Indeed, there are mainly three types 
of dissipation sources for heat engines: friction, heat leaks 
and finiteness of the heat transfert rate between the heat 
reservoirs and the engine. For a heat engine to produce 
power, the working fluid, which transports heat, has to 
move; such displacement of working fluid generates a heat 
flux proportionnal to the matter flux. Hence, any addi- 
tionnal exchange of heat between temperature reservoirs 
can only result in a loss: energy flows through the sys- 
tem without any additional production of power. As dis- 
cussed in Ref. [isl, heat leaks cannot be used to ensure 
causality: they exist independently of the working fluid 
displacement and, as such, they may be considered as 
a parallel parasitic process. For endoreversible models, 
causality is due to finiteness of the heat transfert rate. 
Friction as source of irreversibility and its consequences 
on EMF are discussed in detail in Sec. IVII 

To demonstrate that the CA efficiency is, as a mat- 



ter of fact, an upper bound for EMP, reached only when 
strong coupling is imposed. Van den Broeck used a cas- 
cade construction of heat engines in series [14]. Each 
generator of this chain placed between two heat reser- 
voirs at temperatures Thot and Tcoid respectively, is as- 
sumed to be working at maximum output power; this 
type of setup was further investigated by Jimenez de 
Cisneros and Calvo Hernandez [19, 20]. These previous 
works considered generic heat engines; in the present ar- 
ticle, we focuse on a particular type of system: thermally 
coupled thermoelectric generators. As emphasized by de 
Groot, ''The phenomena of thermoelectricity have always 
served as touchstones for various theories of irreversible 
phenomena ^^ [2l| . 

In a recent work [18], our choice to use a TEG as a 
model system to examine the impact of irreversibilities 
sources on EMP, proved fruitful. We demonstrated that 
the CA efficiency is specific to endoreversible engines and 
hence cannot be used to characterize engines with inter- 
nal dissipation. The fact that our results contradict the 
main message of Van den Broeck in Ref. [14], i.e. the 
universality of the CA efficiency for engines working at 
maximum power in the strong coupling regime, provided 
the impetus to investigate the behavior of a chain of ther- 
moelectric generators in the strong coupling regime in 
order to uncover and explain the subtleties at the heart 
of the physics of energy conversion efficiency of engines 
operating at maximum output power. 

Contrary to the works presented in Refs. [Ij, UH, [20| 
we do not consider an infinite chain but rather a simple 
(though nontrivial, as it turns out) cascade construction 
made of only two heat engines. Our article is organized 
as follows. In Sec. [Ill we use a numerical model to study 
the power and efficiency of the whole system depending 
on the particular electrical load for each TEG. In Sec. lIIH 
using the conce2t of effective thermal conductance intro- 
duced in Ref. [22| , we derive an analytical expression for 
both maximum power and efficiency at maximum power. 
We go further in the analysis of the EMP using our nu- 
merical model in Sec. IIVI where we consider cases when 
the two TEGs characteristics are very different in order 
to study the impact of dissipation localization. Sec. |Vl 
is devoted to internal dissipation. A discussion of our 
findings followed by concluding remarks end the article. 
Further detail on strong flux coupling and Joule heating 
are provided in a series of appendices. 



II. SIMULATION OF TWO TEG THERMALLY 
COUPLED IN SERIES 



The study of a system composed of two thermoelectric 
modules thermally coupled in series is sufficient to pro- 
vide insight into the behavior and main characteristics of 
a cascade construction of coupled heat engines. 
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FIG. 1. Thermoelectric (left) and thermodynamic pictures of 
the thermoelectric generators thermally in series. 



necessary to reach the Carnot efficiency and recover the 
Curzon-Ahlborn value at maximum power for some spe- 
cific working conditions [14]. 

Consider the association of two TEGs thermally in se- 
ries as shown in Fig. [H The quantities defined above, 
now associated with each generator, are denoted with 
subscripts 1 and 2 respectively. One of the main hy- 
potheses we make for this composite model system is the 
adiabaticity of the connection between the two TEGs: 
all the heat released from TEGi is entirely injected in 
TEG2; there is no heat loss at the junction. This relation 
of continuity is the key point to derive the intermediate 
temperature, T^, at the junction point between the two 
TEGs. 

The heat flux rejected by TEGi reads: 



r(out) _ ^ rp T Rih 
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(3) 



Application of Ohm's law yields the electrical current in 
TEGi: h = ai{Thot - Tm)/{Ri + Rn), so we obtain : 



A. Framework 

Before considering the whole system, we give the phe- 
nomenological relations that define the behavior of each 
generators. The individual behavior of each TEG, placed 
between two temperature reservoirs at Thot and Tcoid, can 
be described using a force-flux formalism, relating the 
electrical and thermal currents, / and Iq to the electri- 
cal potential and temperature differences, AV and AT, 
across the generator |22|]: 
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Each TEG is characterized by its isothermal electrical re- 
sistance i?, its thermal conductance at zero electrical cur- 
rent (open circuit )i^^_Q, and its thermopower (or Seebeck 
coefficient) a. The quantity T is the average temperature 
inside the TEG, which we take as T = (Thot + Tcoid)/2. 
The generator performances are characterized by the so- 
called figure of merit ZT defined by ZT = aT/{RK^^^). 
Heat is transported by conduction inside the TEG but 
advection also is a process to account for in presence of 
a nonzero electrical current [22] since transport of heat 
directly is related with the motion of charge carriers in 
the particular direction set by the flow of electric charges 
in the system. Both heat transport processes may be as- 
sociated into an equivalent thermal conductance Kteg 
for the generator that includes advective and conductive 
terms. In the present work we consider TEGs in the 
strong flux coupling regime: the electrical current and 
the thermal current are proportionnal. This condition 
implies that the thermal conductance at zero electrical 
current is zero and that i^xEG reduces to the advective 
part of the thermal conductance (a more exhaustive ex- 
planation is given in Appendix A). This assumption is 
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The additional term associated with Joule heating is dis- 
cussed in the Appendix C. Similarly, the incoming heat 
flux inside TEG2 is given by: 
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Continuity of heat current, Iq^^^ = 1qI\ yields a second- 
order equation in T^ (which we solve numerically for a 
large range of values of the load electric resistances R^i 
and Ri2)' The knowledge of the intermediate temper- 
ature naturally leads to the values of both power and 
efficiency for the whole system as well as for each TEG: 
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The total output power P simply is the sum of Pi and 
P2. The conversion efficiency is defined by : 
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where Iq is given by: 
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FIG. 2. Output powers as functions of the load resistances 
scaled to the internal resistances R^i/Ri and RgilRi'- panel 
a for TEGi; panel b for TEG2; panel c for the whole system. 



B. Results 

The electrical output powers produced and the related 
efficiencies (scaled to the Carnot efficiency, t^c) for each 
thermogenerator and for the whole system as functions of 
the load resistances Rg^ and Rg2^ are plotted in Figs. (J2j) 
and (|3]). For these simulations we took fixed values for 
the temperatures: Thot = 301 K and Tcoid = 299 K; the 
TEGs are similar: Ri = R2 = 0.01^^, and ai =0^2 =0.1 
mV K~^. For the sake of clarity the load resistances are 
normalized to the coresponding TEG internal resistance. 

A first striking result is the appearence of a whole set 
of couples of load resistances {Rii^Ri2) leading to the 
maximum power: each couple corresponds to a particu- 
lar working condition. Another interesting result is that 
the two generators never work at maximum power at the 
same time. We even observe the existence of an interme- 
diate range within which both generators are far from op- 
timal performances. However the maximum power of the 
whole system seems to remain constant even in these in- 
termediate working conditions. Hence, for a chain of heat 
engines operating in the strong coupling regime, there is 
no direct relation between the particular working condi- 
tion of each component and the global working condition; 
this confirms the findings presented in Ref. [20]. 

To analyze the efficiency of the chain we use an expres- 
sion obtained first by Harman [23|. For N heat engines 
thermally coupled in series, he demonstrated that the 
total efficiency r] can be expressed as a function of the 
efficiency rji of each engine composing the chain: 




FIG. 3. Conversion efficiencies (scaled to the Carnot ef- 
ficiency) as functions of the load resistances Rn/Ri and 
Ri2/R2- panel a for TEGi; panel b for TEG2; panel c for 
the whole system. Intermediate temperature Tm on panel d. 



The above expression reduces to: 



r] = m^V2- mm 



(10) 



for two heat engines. 

For each generator the plots showing the efficiency as 
function of working conditions can be simply explained: 
for a specific generator, efficiency, unlike power, depends 
only upon the associated load resistance and not upon 
the other generator working conditions. So, if Ri ^ i?, 
i.e. for an open circuit, the efficiency tends to the Carnot 
value (under the assumption of strong coupling) whereas 
if Ri <C R^ i.e. for a short circuit, the efficiency tends 
to zero. From there, the global efficiency trend can be 
obtained from the Harman relation, Eq. (J9j), which we 
use to also check that, as derived by Van den Broeck 

[14] , the efficiency defined by 77 = 1 — ( ^^^ J reproduces 

itself upon concatenation of heat engines. Equation (J9j) 
also shows that if an engine works with zero efficiency, it 
does not modify the efficiency of the whole system: thus 
a chain composed of two engines with one operating at 
zero efficiency and the other one at the Carnot efficiency, 
works with an efficiency equal to the Carnot value. One 
does not need all the engines of the chain to work at 
maximum efficiency to make the chain work at maximum 
efficiency, as stated in [20,]. 



III. ANALYTIC DERIVATIONS 



TV 



77 = 1 






To gain further insight into the behavior of the coupled 
((^\ TEGs, we derive simple analytical expressions for output 
power and efficiency at maximum power. 



A. Derivation of the intermediate temperature Tm 

The calculations that follow are based on the assump- 
tion of small conversion efficiency, which is easily satisfied 
in the linear regime. Adopting the same approach as in 
Ref. [22j, we obtain: 



^TEG2 + ^TEGi 
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with 
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^TEGs = 






c^i(Tn, + Tcoid)/2 
R2 + Rl2 



(12a) 



(12b) 



By making the reasonnable approximation that T ^ 
{T^ + Tcoid)/2 - (Thot + T^)/2, we end up with the 
following analytical formula for T^: 



To determine the working conditions for each generator 
leading to power maximization we solve these two equa- 
tions: 



dRTi 



and 



dP 
d^ 



(16) 



Both lead to the same relation between the electrical re- 
sistances and the Seebeck coeffients: 



aliRi, - Ri) = al{R2 - Ri2), 



(17) 



which refiects the behavior displayed on Fig. Oc: power 
maximization can, in principle, be reached through an in- 
finite number of combinations of working conditions for 
the two systems. Furthermore, it is in good quantitative 
agreement with the simulation. It is surprising that sat- 
isfaction of electrical impedance matching for both TEGs 
always leads to a power maximization, as though the 
TEGs did not experience their mutual thermal infiuence, 
in that particular case. Moreover, under the approxima- 
tions made, the maximum power reads: 



al/{Ri + Rn)Ti,ot + ay{R2 + Ri2)Tcoid 
ay{Ri^Rn)^ay{R2^Ri2) 



(13) 



Equation (p!3|) is useful to better understand Fig. (J3ld), 
and the asymptotic values of the curves, in particular. 
When both TEGs are shorted (which amounts to hav- 
ing Rii <^ Ri and R£2 "C i?2), the temperature Tm is 
the weigthed average of the reservoir temperatures; each 
weight being the power factor of the related TEG. If only 
one TEG is shorted, T^ becomes close to the associated 
reservoir temperature. For example, when only TEGi is 
shorted, R£i <C Ri and R£2 ^ ^2, the temperature T^ 
tends to Thot- 
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which shows that it is independent of the particular val- 
ues of the load resistances as long as Eq. (pT|) is satisfied. 
Now it is instructive to examines the extreme cases 
when either R£i or R£2 tends to zero: if one of the gen- 
erators (for exemple TEGi) is electrically shorted (i.e. 
Rii <C i^i), it can be viewed as a simple thermal re- 
sistance of conductance iCcontact = ^tegi = o^iT/Ri. 
In that case, and still with the strong coupling assump- 
tion, electrical impedance matching is given by R12/ R2 = 
1 + a^T/(i?2i^contact) (scc Eq. (14) of [23). After sim- 
plifications. 



B. Output power 



Knowledge of the analytical form (p!3|) of the tempera- 
ture Tm permits to rewrite the individual output powers 
of the two thermogenerators given in Eqs. (|6a|) and (|6bp : 
the ouput power for TEGi now is given by: 



Ri2 = R2 -\ 2 ^1 



(19) 



and we see that we recover exactly Eq. (p!7|) . The same 
conclusion holds if TEG2 is shorted instead of TEGi; 
only the subscripts need to be swapped. 



Pi = ai(Thot - Pcold) 
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The expression for P2 is obtained simply by swapping the 
subscripts 1 and 2; so the total output power dehvered 
by the two coupled TEGs is given by: 



P = alaliTi,ot-Tcoidy 
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C. Efficiency at maximum poAver 

Let us now turn to the efficiency at maximum power 
^(Tmax), which is given by the ratio of the maximum 



power Pmax to the incoming heat fiux I^ 
pends on P^ax: 
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In the above expression, we see that the dependence of 
the incoming heat flux Iq on Pmax is through the tem- 
perature Tm- To determine Tm(Pmax), we use the ana- 
lytic form of Tm given by Eq. (p!3|) , and the condition 
for electrical impedance matching, Eq. (fT7|) . Then a 
straightforward calculation leads to the following simple 
expression of 7^(Pmax)- 



^(Pmax) = 



m 



(22) 



We thus recover a well known result, except that this is 
not through an expansion of /^(Pmax)- In this section we 
have simplified the problem to fully grasp the subtleties 
of efficiency at maximum power; in the next section we 
base our analysis on the results of the simulation shown 
above in Sec. (|TT|) . 



IV. ANALYSIS OF THE EFFICIENCY AT 
MAXIMUM POWER 

To study the efficiency at maximum power we ffrst 
have to determine the various working conditions corre- 
sponding to this particular engine operation mode. This 
amounts to locating the coordinates of the maxima of 
the curves displayed on Fig. [21 We consider the maxima 
along the lines determined by a constant ratio of R12/R2 
to Rii/Ri. Variation of this ratio yields a description of 
the whole set of possible combinations of working condi- 
tions leading to power maximization [i.e. the maxima of 
Fig. ©]. 

The efficiency associated with the maximum power 
condition is shown on Fig. |4]as a function of the ratio of 
^^2/^2 to Rii/Ri^ i.e. for various combinations of ther- 
moelectric generators electrical load conditions. We con- 
sider three cases: i/ the TEGs are identical; ii/ the inter- 
nal resistance of TEGi is much larger than that of TEG2; 
in/ the internal resistance of TEGi is much smaller than 
that of TEG2. In all cases, the Seebeck coefficients of 
both generators are equal. 

As expected (see e.g. [18, z5]) /^(Pmax) is a little larger 
than 770/2, and thus different. We find an additional 
dependence of this efficiency at maximum power on the 
ratio of internal resistances R2/ Ri. We can make a gen- 
eral statement: a stronger dissipation on the hotter side 
than on the colder side yields a higher efficiency at max- 
imum power. This result is consistent with the case of a 
single TEG with dissipative thermal contacts [18|. 

The condition {R12/R2) ^ (Rii/Ri) means that the 
first TEG is closer to the short-circuit state than the 
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FIG. 4. Reduced efficiency at maximum power depending on 
the ratio of Ri2 and Rn normalized by R1/R2. The three 
curves correspond to the cases where Ri = R2, Ri = 10000 x 
R2 and R2 = 10000 x Ri. 



second one and hence that TEGi dissipates more power 
than TEG2. Furthermore, if Ri ^ i?2, the capacity of 
TEGi to internaly dissipate power is greater than that 
of TEG2, which results in an increase of /^(Pmax) • the 
curve of /^(Pmax) as a function of {R12 / R2) / {Rii / Ri) is 
shifted to higher values. 

We note that for the symmetric case {Ri = i?2), the 
EMP 7^ (Pmax) varies around the value found for the case 
of a single TEG with dissipative contacts to the temper- 
ature reservoirs: 770/2 + Vc/^- ^^^ factor 1/8 at the 
second order of the expansion in r]c of 77 (Pmax) is con- 
sidered as universal for symmetric low-dissipation heat 
engines opera ting in the strong coupling regime: first in- 
tuited by Tu [2J1, this result was demonstrated later by 
Esposito and co-workers [25]. 

When the ratio R1/R2 is changed, 7/ (Pmax) may vary 
between 770/2 + 7^^/16 and 7^0 /2 + 37^^/16 depending on 
the ratio of dissipation between hot and cold sides. This 
result brings back to mind the expression for 7/ (Pmax) 
first derived by Schmiedl and Seifert for a stochastic heat 
engine [26] : 



^Ss(Pmax) = 
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where the parameter 7 characterizes the distribution of 
dissipation. The interest of considering an association 
of two coupled TEGs rather than a single one, as stud- 
ied in Ref. [18], is that 7 may differ from 1/2, a value 
imposed in the case of single TEG by the symmetry of 
the internal disspation mechanism (Joule effect). If each 
TEG still obeys to Joule effect symmetry, a dissymetry 
of either the TEGs themselves {Ri 7^ i?2), or the working 
conditions of each module (through R£i and R£2) leads 
to a significant modification of 7. Combination of TEGs 
seems a simple way to achieve tuning of the parameter 
7- 
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FIG. 5. Heat flux from tfie fiot reservoir (left) and maximum 
power (right) as functions of the ratio of Ri2 to Rn normal- 
ized by Ri/R2- In this case, Ri — R2- 



One may wonder about the nature of the dependence 
of 7?(Pmax) ou the ratio {Ri2/R2)/{Rii/Ri)- If the 
maximum power value is actually constant as derived 
in Sec. piip . the increase in efficiency for {R12/R2) ^ 
{Rill Ri) should be the consequence of the decrease of 
incoming thermal current in the first TEG. To check 
this hypothesis we computed both the maximum power 
value and the associated incoming thermal flux for the 
case where Ri = R2] the results are plotted in Fig. (|5]). 
We see that the maximum power and the thermal cur- 
rent actually depend on the specific working conditions 
of both engines, and that they follow the same trend 
as 7?(Pmax)- This discrepancy with the analytical ap- 
proach where the maximum power is constant may be 
explained by the fact that in the analytical derivation we 
made an approximation for the average temperatures: 
{T^ + Tcoid)/2 « (Thot + Tm)/2 ^ T; this led to a simpli- 
fication of T in Eq. (pT|) . We believe that the small vari- 
ations are, in fact, linked to the ratio of (Tm + Tcoid)/2 to 
(Phot +Pm)/2, which is close to one but varies sufficiently 
with the load resistances to change the conditions for 
power maximization. To be rigorous, we should say that 
the true power maximization is reduced to one working 
condition that the first TEG is shorted, but we consid- 
ered that the maximum power condition is satisfied for 
all local maxima. 

We end this section with one surprising fact arising 
from our analysis of efficiency at maximum power: inde- 
pendently of the specific values of TEGs' parameters (ai, 
0^2, Ri and P2), whenever electrical impedance matching 
is achieved for both TEGs simultaneously {Rn = Ri and 
R£2 = R2) we obtain the value of the Curzon-Ahlborn 
efficiency at the second order in the Carnot efficiency: 
rjc/2 + T^yS, for the EMP of the whole system. For 
this particular working condition, the two TEGs have re- 
spectively the CA efficiency too, which is coherent with 
Eq. (|9|). 



V. ON THE IMPORTANCE OF INTERNAL 
DISSIPATION 

Energy dissipation is necessary to ensure the causal- 
ity of the processes during the operation of heat engines 
flSj . For endoreversible engines as the one presented in 
Ref. [10], each flux is related to generalized forces through 
linear relations, in the engine as well as in the contacts 
to the reservoirs. Inclusion of internal dissipation in the 
model obviously breaks the assumption of endoreversibil- 
ity and necessary implies addition of the dissipated en- 
ergy in the heat balance. Since the mathematical ex- 
pression for the quantity of dissipated heat is given by 
the product of a generalized force by a flux, the addi- 
tional term corresponding to internal irreversibilities is 
a quadratic term simply because of the proportionnal- 
ity between force and flux. Therefore, even though the 
analysis of the exoreversible engines may have to be per- 
formed within the framework of linear irreversible ther- 
modynamics, the quadratic term cannot be neglected 
since it is necessary to ensure the causality of the pro- 
cesses associated to the operation of the engines. There- 
fore, though the dissipative term is not stricto sensu a 
linear form, we propose to consider it as part of the lin- 
ear irreversible thermodynamics nonetheless, as a direct 
consequence of the hnearity of the constitutive equations 
of the model. A similar viewpoint was already retained 
by C alien and Welton in Ref. [23 • Recent publications 
[28|, [29] addressed this issue of linearity in irreversible 
thermodynamics. In particular Izumida and Okuda laid 
the foundations of an extended Onsager model consider- 
ing the supplementary internal dissipative term |28j. 

The following discussion on the importance of the in- 
ternal dissipation provides further insight into the impact 
of dissipation localization on EMP. In Ref. [18] we saw, 
while studying dissipative thermal contacts, that dissipa- 
tion occuring on the hot side of the engine enhances its 
EMP. The analysis of Fig. (HI shows that the same trend 
is recovered here for two thermally coupled TEGs: as 
already discussed in Sec. (jIVj) a dissipation located pref- 
erentially on the hot side, with for example Ri > R2, al- 
lows the system's EMP to overcome the Curzon-Ahlborn 
efficiency bound (or more precisely, the Schmiedl-Seifert 
efficiency in that case). Note that the idea of heat inter- 
nally dissipated running back to the hot reservoir, hence 
allowing better performance of the heat engine was re- 
cently discussed in a paper on the thermodynamics of 
hurricanes [30]. The author stated that these meteoro- 
logic phenomena possess the capacity to preferentially 
dissipate energy in the higher temperature areas, which 
in turn enhances their efficiency. Though we will not 
comment on the derivations made to characterize the hur- 
ricanes, this idea of efficiency enhancement due to energy 
dissipation in the hotter region of a thermodynamic sys- 
tem is consistent with our findings on the behavior of a 
chain of TEGs. 

To extend our discussion we consider an irreversible 
Carnot engine. As this engine experiences friction when 
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FIG. 6. Temperature vs. Entropy for different engine models 



tfie piston moves, energy is dissipated as fieat is pro- 
duced. However, as tfie cfiamber is supposed to be adia- 
batic, fieat remains confined inside and may only transit 
to tfie reservoirs during the isothermal steps of the cycle. 
The heat produced during the adiabatic step from Tcoid 
to Thot may thus be used during the next cycle. This is 
analogous to the Joule heating flowing back to the hot 
reservoir, thereby reducing the incoming heat flow. Adi- 
abatic expansion actually is the key point to understand 
the importance of the internal disspation inside an exore- 
versible Carnot engine. 

A temperature-entropy (T — S) diagram highlights the 
fact that the adiabatic expansions are not always isen- 
trop ic for a non-ideal Carnot engines as shown in Ref. 
[3IJ . The behavior of three different configurations for 
the non- ideal Carnot engine is depicted on Fig. ([6|): on 
the left side is the exoreversible case where all dissipa- 
tion takes place inside the engine; in the middle, is the 
endoreversible case studied by Novikov, and Curzon and 
Ahlborn; on the right side is a mixed situation involving 
irreversibilities that occur both internally and externally. 
The endoreversible cycle is the easiest one of the three to 
understand: only the temperature of the isothermal steps 
are modified due to the presence of dissipative thermal 
contacts. For the exoreversible engine, isothermal steps 
occur at Thot and Tcoid respectively; the effect of non- 
idealities lies in the adiabatic steps: as frictions produce 
heat, the entropy increases during these steps, but be- 
cause of the hypothesis of adiabaticity this heat remains 
in the chamber, the main consequence being that the ad- 
ditionnal heat has to be released into the heat reservoirs 
during the isothermal steps. This additionnal heat is rep- 
resented on Fig. (|6|) by the purple triangular area for the 
step from cold to hot, and by the blue triangular area the 
step from hot to cold. However these two quantities have 
totally different impact on the engine performances: the 
heat produced during the cold to hot step contributes to 
reduce the absorbed heat during the isothermal step at 
Thot and hence allows to increase the efficiency. On the 
contrary, the heat produced during the hot to cold step 
is released at the cold side and cannot be recycled: it is 
only loss. 

Now, considering both internal and external dissipa- 
tions, we account both for effective temperature relax- 
ation and heat production during the adiabatic steps. In 
Ref. [31], Wu and Kiang considered this non trivial case. 
They introduced a quantity R to take account of the ratio 



of useful heat production in cold to hot step and of pure 
loss heat production in hot to cold step. This quantity 
possesses similarities in its construction with the param- 
eter 7 appearing in the formula derived by Schmiedl and 
Seifert, Eg. (|23|) . Unfortunately the ratio R as defined in 
Ref. [31| turns out to be innacurate as Wu and Kiang ne- 
glected the temperature modification at the engine edges: 
they considered that the temperature of the isothermal 
steps remains the same as that obtained for an endore- 
versible engine, but the heat produced during adiabatic 
steps actually contributes to a slight increase of these 
temperatures when frictions are included. This tempera- 
ture relaxation is similar to the mean temperature mod- 
ification of a TEG when shifting from endoreversible to 
exoreversible conditions as illustrated in [18]. 

Despite its crucial importance, the adiabatic step often 
is neglected in the analysis of heat engines. Furthermore, 
we want to emphasize the fact that Novikov had already 
considered the fact that internal dissipation plays an im- 
portant role in the behavior of heat engines. In Ref. [3], 
he introduced a parameter A (quite similar to the param- 
eter 7 of Schmiedl and Seifert, and the ratio R in Wu's 
and Kiang's work) which characterizes the distribution of 
the entropy produced by internal irreversibilities between 
the two adiabatic steps. However in his article, Novikov 
chose to consider only an endoreversible model and con- 
sequently neglected the influence of the parameter A. 

The EMP bound as defined by Gaveau and co-workers 
[32] is reached for the extremal values of R: in the less 
favorable case, all irreversibilities occur during the hot- 
to-cold stage, and the EMP is only r]~ = 7^c/2; whereas 
in the most favorable case, when all irreversibilities oc- 
cur during the cold-to-hot stage, the EMP may reach 
77+ = 77c/ (2 — 77c). In the case of a single TEG, as heat 
is produced by the Joule effect, it is equally distributed 
between hot and cold heat reservoirs, thus leaving no de- 
gree of freedom to tune significantly the EMP: one can 
only shift from endoreversible to exoreversible conditions. 
Nevertheless, for a tandem setup as the one considered 
in the present article, we can break this natural sym- 
metry by changing both individual working conditions 
and TEGs parameters. These tuning parameters indeed 
lead to EMP modifications as shown in Fig. (j4]). How- 
ever, the bounds r]~ and 77+ cannot be reached. Since 
a thermogenerator must produce Joule heating in order 
to deliver power, it is impossible to fully release the pro- 
duced heat to a single heat reservoir; this explains why, 
in our case, the EMP is comprised between 7/0/2 + 77^/16 
and r]c/2 -h 37/^16. 



VI. DISCUSSION AND CONCLUSIVE 
REMARKS 

Our viewpoint For a tandem construction such as the 
one proposed in Ref. [IJ], there exist several options to 
tune the working condition of the whole system for given 
heat reservoir temperatures. Indeed the junction tem- 



perature T^ and the working conditions for each en- 
gine composing the tandem, form three different lever- 
ages that can be used for such purpose. However, one 
should keep in mind that these cannot be employed all 
at the same time: internal laws governing the system, 
e.g. heat flux continuity at the interface between en- 
gines, determine the third parameter when the other two 
are fixed. For pratical applications, imposing the elec- 
trical load resistances seems the simpler one. Therefore, 
our point of view is that the temperature profile should 
not be fixed, but rather deduced from the working con- 
ditions of each component of the chain; this view differs 
from those given in Refs. [1^, [19| and [20]. Interestingly, 
we recover the fact that, as shown by Van den Broeck 
\\A\ , when both TEGs work with the Curzon- Ahlborn ef- 
ficiency, the whole system is working at maximum power 
with Curzon- Ahlborn efficiency too. However, we have 
demonstrated that in that case, surprisingly, none of the 
TEGs works at maximum power. 



Difference with compatibility approach In the com- 
patibility approach of thermoelectricity [3J-[36|], the 
problem is quite different as the TEGs are not only ther- 
mally in series but are electrically coupled too. Thus, 
one cannot tune the working conditions of the TEGs in- 
dependently: we loose a degree of freedom as only one 
load resistance, common to the two TEGs, remains, as it 
is the case in Ref. [23|. 



Conclusion Using a particular tandem construction 
made of two thermoelectric generators thermally in series 
but electrically independent, we have derived general re- 
sults concerning heat engines with a special focus on the 
efficiency at maximum power. First, we have obtained 
the working conditions yielding power maximization of 
the whole system: we noted that this condition is not 
unique and that it can be obtained for various combina- 
tions of the electrical load resistances of each generator. 
Through analytical derivations, we obtained an expres- 
sion that defines the couples of resistances maximizing 
the power. Yet, the analytical approach was not sufficient 
to fully describe the variations of efficiency at maximum 
power; so we turned to numerical simulations. We recov- 
ered the fact that the Curzon Ahlborn efficiency can be 
overcome for some specific cases. We then explained the 
variations of EMP focusing on the importance of inter- 
nal dissipation. We demonstrated that the distribution 
of this internal dissipation is the key point to understand 
issues related to the EMP: if this dissipated heat flows 
back to the hot heat reservoir, it can be viewed as a re- 
cycling process, i.e. one has not to supply this heat dur- 
ing the next isothermal process, whereas dissipated heat 
that flows to the cold heat reservoir, amounts to pure 
loss. The tandem construction allowed us to vary this 
distribution and demonstrate its impact on the efficiency 
at maximum power. 
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Appendix A: Thermoelectric module operating in 
the strong flux coupling regime 

Adopting a local description of the operation of a 
thermoelectric module, the coupled transport of electric 
charges and heat may be expressed as follows in the force- 
flux formalism: 
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where Lij are the kinetic coefficients associated with the 
phenomenon. 

It is often convenient to rewrite the fluxes as a function 
of the gradients of temperature T and electrochemical 
potential /i: 
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where a is the Seebeck coefficient, a the electrical con- 
ductivity, i<ij._^ the thermal conductivity under open cir- 
cuit condition and e the elementary electric charge. To 
evaluate the coupling strength between the electrical and 
thermal fluxes we introduce the parameter q defined as: 



Li2 



yLiiL 



(A3) 



22 



After substitution of the kinetic coefficients by their ex- 
pressions in the particular case of thermoelectric trans- 
port (see Ref. [37]) we obtain: 
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where sgn denotes the signum function. At first sight, 
a large value of the figure of merit ZT is sufficient to 
ensure a strong flux coupling inside the TEG. However, 
we should keep in mind that the strong coupling also is 
characterized by a relation of proportionnality between 
electrical and thermal fluxes: even the limit ZT — )► oo 
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FIG. 7. Power (first column) and efficiency (second column) 
for different values of |^| (or equivalent ly ZT). a-b: |^| = 1 
(ZT = oo), c-d: 1^1 = 0.87 {ZT = 3), e-f: |^| = 0.05 {ZT = 
0.003). 



is insufficient on its own. The strong coupling regime 
may be obtained with /^/^^ = 0; in other cases there still 
is an additionnal term in the expression of the thermal 
current: Jq = aTJ — nj^^VT. Furthermore, assuming 
a ^ oo necessarily implies that J = since it is the only 
way to avoid a divergence of the thermal flux. 



Appendix B: Case with no strong flux coupling 
assumption 

To analyze the influence of the coupling parameter q 
between the electrical and thermal flux, we choose to tune 
it keeping constant internal resistance and Seebeck coef- 
ficient and only changing the value of the thermal con- 
ductance under open-circuit condition K^^^. For sake of 
simplicity, the two generator are always identical. The 
impact of the q parameter is evidenced in Fig. ([7|). For 



strong coupling, i.e. g = 1, the influence of each load re- 
sistance is strongly correlated to the working condition of 
the other TEG. When this hypothesis is relaxed and ter- 
mal leaks are raised, i.e. K^^^ increased, the two TEGs 
working conditions become more and more independent. 
When the advective part of the thermal flux, controled by 
electrical load trough the value of the current, becomes 
negligeable compared to the conductive heat flow the per- 
formance of the whole system are given only by the ad- 
dition of the two TEGs' taken independently: there is no 
more correlation between the two subsystems as shown 
in the third line of Fig. ([7]). 



Appendix C: Comments on Joule heating 

To satisfy the steady-state heat equation, sometimes 
referred to Domenicali's equation [17, 38], one needs to 
add in an ad hoc fashion the heat generated by the Joule 
effect to the thermal flux balance. Solving v ■ Jq = 
( V /i) • J, we obtain: 
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Assuming constant system parameters {va = and 

V K,^_^ = 0), and a constant temperature gradient (given 
by the ratio of the temperature difference to the length of 
the module) as in the Onsager-Callen formalism |15l, fig, 
the left hand side of Eq. ()Cip reduces to zero. This yields 
J'^ /a = 0, which is impossible in the presence of a finite 
electrical current in the system. To satisfy the heat equa- 
tion, an extra term must be added to the heat flux. 

The force-flux formalism only deals with average quan- 
tities inside the heat engine. To account for the differ- 
ence between input and output fluxes the thermal flux 
may be made asymmetric by addition of algebraic terms 
characterizing the total Joule heating, in the incoming 
and outgoing fluxes. Each end of the system is allo- 
cated half the heat as demonstrated in Ref. |39|. This 
equipartition remains valid at the mesoscopic scale |4Q| |. 
One way to avoid the ad hoc addition of Joule heating 
is to allow the temperature gradients on the system's in- 
put and output sides to be different: to do so one may 
define an internal temperature. By solving numerically 
both Onsager-Callen relations and the heat equation, it 
is easy to show that one does not actually need addition 
of Joule terms: the system's internal temperature reflects 
the Joule heating [41.] . 
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